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FOREWORD 


This  report  examines  a least  squares  method  for  the  solution  of 
orbital  and  bias  parameters  for  multiple  Global  Positioning  System 
satellites.  The  method  includes  the  solution  for  timing  errors  present 
in  the  system. 
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ABSTRACT 


A least  squares  method  is  presented  for  the  solution  of  orbif-1  and 
bias  parameters  for  multiple  Global  Positioning  System  satellites  being 
observed  by  a common  station  array.  The  effects  of  radiation  pressure 
modelling  errors  on  timing  solutions  are  examined  for  both  multiple  and 
single  satellite  processing  using  a time  equation  consisting  of  a bias 
and  linear  drift  term.  Results  indicate  that  the  errors  in  orbital  and 
timing  solutions  caused  by  modeling  errors  in  radiation  pressure  are 
reduced  by  the  multiple  satellite  processing  scheme.  However,  adequate 
time  solutions  from  single  satellite  processing  are  achievable  if  accu- 
rate radiation  pressure  modeling  is  available. 
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i . n.TKonrcTioN 


r>-e  N'.ivstar  Global  Positioning  System  (CPS)  is  a passive  all  weather 
navigation  satellite  system  proposed  for  the  1980-1990  time  frame.  The 
system  is  currently  in  the  prototype  evaluation  phase  with  two  experi- 
mental models  to  be  launched  in  the  next  three  years.  These  prototype 
models  are  being  designed  and  built  by  the  Naval  Research  Laboratory  to 
test  further  the  concept  of  passive  satellite  navigation  based  on  highly 
accurate  frequency  standards.  Timation  IH-A  is  the  first  in  this  series. 

'..•hen  operational  the  CPS  system  will  enable  users  to  determine  their 
three  dimensional  position  and  time  instantaneously.  Ranging  measure- 
ments taken  simultaneously  from  at  least  four  satellites  will  be  reduced 
to  determine  these  parameters.  Anticipated  positional  accuracies  in  the 
horizontal  and  vertical  axes  are  better  than  10  meters  90  of  the  time. 

obviously  there  will  exist  a critical  need  for  precise  orbit  deter- 
mination and  prediction,  and  for  extremely  stable  oscillators  aboard  the 
satellites.  The  CPS  system  as  envisioned  will  employ  atomic  frequency 
standards  In  each  satellite  to  assure  high  oscillator  stability.  The 
'■aval  Research  Laboratory  is  currently  developing  the  clock  technology 
necessary  to  -*et  CPS  accuracy  goals.  Prototype  models  (Timation  series) 

will  test  the  feasibility  of  using  rubidium,  then  cesium  standards  in  a 
spicv  environment. 

Current  plans  call  for  approximately  five  operational  ground  stations 
to  track  all  system  satellites  for  orbit  determination.  These  tracking 
stations  will  also  be  equipped  with  high  stability  oscillators.  But  even 
with  atomic  standards  there  are  offsets  and  drifts  in  frequency. 
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Since  frequency  (timing)  errors  at  the  comoon  tracking  stations  will  be 
reflected  in  the  computed  orbits  of  all  GPS  satellites  a procedure  can 
be  developed  which  will  use  datn  from  vaiious  satellites  to  determine 
system  ephemerides  and  timing  errors  precisely.  Such  a procedure  is 
examined  in  this  report.  The  method  models  the  timing  error  of  each 
system  clock  as  a bias  and  linear  drift  term.  The  assumption  is  that 
this  approximation  will  be  accurate  for  atomic  standards  over  a few  day 
period.  This  assumption  is  supported  by  rubidium  frequency  data  taken 
at  the  Naval  Research  Laboratory  (NRL)  r Reference  1 . The  method  may 
incorporate  more  sophisticated  clock  modeling  as  warranted. 

However,  precise  system  time  determination  depends  on  the  accuracy 
of  the  dynamical  model  used  in  orbit  determination.  It  has  been  demon- 
strated that  radiation  pressure  modeling  for  GPS  altitude  satellites  is 
critical  Reference  2\  A radiation  pressure  bias  of  one  percent  may 
introduce  orbit  errors  of  four  meters  during  two  day  orbit  predictions. 
The  effects  that  radiation  errors  may  have  on  the  timing  solutions  for 
CPS  satellites  are  also  considered  in  this  report. 


!!•  ML*1.TI  PLE  SATELLITE  PROCESSING 


!!ie  concept  of  nultiple  satellite  processing  for  the  GPS  system  is 
based  on  the  fact  that  timing  errors  at  common  tracking  stations  will 
affect  the  computed  orbits  of  all  satellites.  Theoretically  a multiple 
satellite  processing  scheme  should  predict  the  parameters  of  the  time 
equation  better  since  it  utilizes  all  observations.  These  timing 
parameters,  which  are  modeled  here  as  a bias  and  linear  drift,  replace 
range  and  range  drift  pass  parameters  associated  with  range  data 
(Appendix  B1 . Knowledge  of  timing  errors  at  one  station  is  required  to 
make  timing  parameters  linearly  independe  ;.  A better  knowledge  of 

system  timing  errors  w:il  allow  a better  determination  of  satellite 
orbital  constants. 

A mathematical  description  of  least  squares  processing  of  multiple 
satellite  data  is  described  in  Appendix  B.  The  procedure  involves 
forming  the  normal  equations  for  each  pass  of  satellite  data,  formally 
eliminating  the  pass  dependent  parameters,  combining  normal  equations 
and  solving  for  orbital  and  liming  corrections. 
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ITT.  S KNSITIVITY  ANALYSIS 


; A.  Radiation  Pressure  Modeling 

j - “ 

An  analysis  was  conducted  to  determine  what  effects  errors  in 
radiation  pressure  modeling  would  have  on  timing  solutions  produced  by 
both  individual  and  multiple  satellite  processing.  The  results  for  all 
cases  assume  the  following: 

(i)  System  clocks  are  accurately  modelled  with  a bias  and 
linear  drift. 

(ii)  Tracking  station  positions  are  perfectly  known. 

(iii)  Clock  errors  are  perfectly  known  for  one  station. 

1.  An  Cnmodelled  Bias 

a.  Multiple  Satellite  Processing 

Pesults  have  been  obtained  using  synthetic  range  data 
from  six  satellites  (Table  1)  and  five  tracking  stations  (Table  2).  The 
standard  error  on  the  range  data  was  50  centimeters.  Again  the  timing 
errors  at  one  station  were  assumed  known  to  a high  degree  of  accuracy  in 
order  to  make  the  timing  parameters  independent.  Cases  were  considered 
with  a radiation  pressure  bias  of  one  and  five  percent.  No  attempt  was 
made  to  solve  for  this  error. 

For  each  case  orbital  and  timing  solutions  were  made 
over  a two  day  period  using  the  least  squares  approach  discussed  above. 
The  results  obtained  for  the  time  bias  and  drift  terms  are  given  in 
Tables  3 and  4 respectively.  The  average  standard  e *ror  on  the  timing 
solutions  is  also  given.  Table  5 gives  the  maximum  difference 
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between  the  true  clock  time  and  the  time  predicted  using  the  multiple 
satellite  solution.  This  difference  is  approximately  a linear  function 
of  radiation  modeling  error  magnitude  when  the  error  appears  as  a bias. 


SATELLITE  ORBITAL  ELEMENTS 
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TRACKING  STATION  NET 


Station 

Number 

Latitude 

Longitude 

Panama 

1 

9.0° 

280.0° 

Florida 

2 

25.6° 

279.6° 

Maryland 

3 

38.7° 

283.5° 

Seychelles 

4 

*4.7° 

55.5° 

Samoa 

5 

-14.3° 

189.3° 

TABLE  2 


« 
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TIME  BIAS  SOLUTION 
vs. 

RADIATION  MODEL  ERROR 

Radiation  Bias 


Satellite 

True  Bias 

n 

5% 

1 

10.0* 

10.0* 

10.2* 

2 

15.0 

15.0 

14.9 

3 

-5.0 

-4.5 

-2.6 

4 

20.0 

20.5 

22.5 

5 

25.0 

25.2 

26.0 

6 

-20.0 

-19.8 

-19.1 

Station 

1 

5.0 

4.9 

4.5 

2 

10.0 

9.9 

9.3 

3 

-10.0 

-10.2 

-11.1 

4 

15.0 

14.6 

12.8 

5 

0.0** 

0.0 

0.0 

Standard  £ 

*ror^l  Sigma>  On  Time 

Bias  Solution:  .07 

nanoseconds 

* nanoseconds 

—station  5 parameters  constrained 


TABLE  3 
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TIME  DRIFT  SOLUTION 


vs , 

RADIATION  MODEL  error 

Radiation  Bias 


Satellite 

True  Drift 

n 

57c 

1 

40.0* 

40.0 

40.3* 

2 

25.0 

25.2 

26.2 

3 

30.0 

29.7 

28.4 

4 

-25.0 

-25.2 

-26.1 

5 

30.0 

30.0 

29.8 

6 

-20.0 

-20.0 

-20.4 

Station 

1 

25.0 

25.1 

25.6 

2 

30.0 

30.0 

29.7 

3 

35.0 

35.0 

35.0 

4 

-20.0 

-19.8 

-19.0 

5 

0 . O’- * 

0.0 

0.0 

Standard  Errorfl  Sigma)  On  Time  Drift  Solution: 

.03  nanosec/day 

* nanoseconds /day 
’^station  5 parameters 

constrained 

TABLE  4 
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MXIMl’M  TIME  SOLUTION  fppop 
M'E  TO  RADIATION  BIAS 


F 


i 


During  Fit  Span 

Radiation  Error  i*, 


A t 

5 Z 

Satel li te 

1 

o.o-’* 

.8* 

9 

.4 

2.3 

3 

.5 

2.4 

<4 

.5 

2.3 

5 

.2 

1.0 

b 

.2 

.9 

St nt ion 

1 

.1 

.7 

2 

.1 

1.3 

3 

.2 

1.1 

4 

.4 

2.2 

0.0 

0.0 

nanoseconds 


TABLE  5 
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b.  Individual  Satellite  Processing 


I’sing  the  synthetic  range  data  generated  for  the 
multiple  satellite  case  individual  satellite  orbit  solutions  were  deter- 
mined which  best  fit  the  data  in  the  least  squares  sense.  Cases  were 
again  considered  where  radiation  biases  of  one  and  five  percent  were 
present.  In  this  processing  mode  timing  errors  were  solved  for  as  pass 
bias  parameters.  The  time  bias  present  in  the  satellite  and  station 
clock  were  combined  into  a range  bias.  Time  drifts  were  likewise  con- 
sidered as  a range  drift.  Again  unmodelled  radiation  bias  resulted  in 
errors  in  the  range  bias  solution  which  were  linear  with  radiation  error 
magnitude.  Table  6 gives  the  average  of  the  pass  timing  errors  due  to 
radiation  bias.  This  value  in  nanoseconds  is  equivalent  to  the  error  in 
computed  range  bias  due  to  radiation  error. 

c.  Comparison  of  Processing  Modes 

Since  each  processing  method  treats  the  timing  para- 
meters differently  a direct  comparison  of  results  is  difficult.  How- 
ever, a comparison  of  Table  5 and  6 can  reveal  information  about  the 
size  of  errors  due  to  radiation  bias.  For  individual  satellite  pro- 
cessing the  average  time  error  for  a pass  is  about  1.2  nanoseconds  for 
1'  radiation  bias  and  6.6  nanoseconds  for  5 7 bias  (Tabl  6).  These 
error  levels  are  for  a pair  of  clocks.  However  in  m»i1wiple  processing 
the  maximum  error  levels  for  the  same  radiation  biases  are  smaller.  If 
the  worst  satellite  and  station  errors  in  Table  5 are  summed  the  result 
is  waller  than  individual  processing  results. 
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INDIVIDTAL  PROCESSING  MODE 
AVERAGE  TIME  ERROR  PER  PASS 


Radiation  Bias 

11 

52 

Average  Pass  Time  Error 

1.2 

6 . 6 

( nanoseconds  > 


TABLE  6 
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Another  way  to  check  on  the  relative  accuracy  of  the  two 
methods  is  to  compare  the  orbit  solutions.  Figure  2 graphically 
illustrates  various  trajectory  differences  resulting  from  a 5 radiation 
bias.  Table  7 lists  each  case  considered. 

In  cases  t and  II  both  radiation  bias  and  timing  errors 
were  present  in  the  data.  Case  III  however  has  no  timing  errors  pre- 
sent. In  no  case  was  radiation  a parameter  of  fit. 

From  Figure  2 it  can  be  seen  that  with  a bias  in  radiatbn 
pressure  modeling  the  multiple  satellite  processing  mode  computes  an 
ephemeris  which  is  more  consistant  with  that  produced  from  data  with  the 
same  radiation  error  but  without  timing  errors.  The  coupling  between 
ttie  range  bias  (timing  errors')  and  the  induced  error  due  to  radiation 
bias  causes  a "poorer”  determined  orbit  with  the  individual  processing 
mode. 

2.  An  Approximate  Radiation  Pressure  Model 

Since  radiaLion  pressure  modeling  is  critical  for  CPS 
altitude  satellites,  attempts  are  being  made  to  model  radiation  pressure 
to  better  than  1 . At  the  Naval  Sur  face  WeaponsCenter  various  radiation 
pressure  models  are  being  developed  for  the  Timation  1 1 1 -A  satellite. 

Hie  asymmetry  of  this  satellite  prohibits  a spherical  model  with  fixed 
surface  to  mas^  ratio  to  be  used  in  orbit  determination.  A model  incor- 
porating varying  exposed  surface  area  and  approximating  the  reflection 
and  absorbtion  properties  of  the  satellite  is  required. 

To  determine  the  errors  that  will  be  introduced  into 
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timing  solutions  due  to  an  approximate  surface  model,  an  orbit  was  gen- 
erated using  a rectangular  box  model  with  sides  of  varying  surface  area, 
reflectivity  aad  absorbtion.  A least  squares  fit  was  made  using  range 
data  with  timing  errors  to  determine  the  satellite's  ephemeris  and  pass 
timing  errors.  The  radiation  model  employed  in  the  orbit  solution  was  a 
sphere.  A scaling  parameter  was  determined  for  this  model.  The  average 
pass  timing  error  introduced  by  using  the  spherical  surface  model  was  .6 
nanoseconds . 

Thus  it  is  likely  that  a radiation  pressure 
nodel  which  approximates  the  physics  of  the  satellite's  orbit  in  terms  of 
radiation  pressure  may  yield  acceptable  timing  results  with  individual 
satellite  processing.  Incorporating  this  model  into  the  multiple 

satellite  processing  scheme  should  yield  very  acceptable  results. 

B . Station  Position  tncertainty 

The  preceeding  results  were  based  on  the  assumption  that 
station  position  locations  are  perfectly  known.  Increasing  the  uncer- 
tainty in  station  position  directly  effects  the  standard  errors  on  the 
time  bias  and  drift  solutions.  Table  8 demonstrates  increasing  standard 
errors  of  the  timing  solutions  as  a function  of  increasing  uncertainty  in 
station  position  for  all  five  stations  considered  in  the  multiple  satel- 
lite case.  Tliis  correlation  between  the  uncertainty  in  station  position 
and  standard  errors  on  timing  solutions  necessitates  accurate 
determination  of  tracking  locations  tor  the  HPS  system. 
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TIMING  SOLUTION  UNCERTAINTY 


vs. 


STATION  POSITION  UNCERTAINTY 


Position  Uncertainty 
(meters) 


Bias  Uncertainty 
(nanosec) 


.07 


.4 


1.0 


TABLE  8 
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Drift  Uncertainty 
(nanosec /day) 


.03 


.3 


.7 


TV.  CONCLUSIONS 


This  report  has  examined  a least  squares  method  for  the  solution  of 
orbital  and  bias  parameters  for  multiple  Global  Positioning  System  satel- 
lites being  observed  by  a common  station  net.  The  method  includes  the 
solution  tor  time  biasing  parameters  present  in  the  system.  The  effects 
of  radiation  pressure  modeling  errors  and  station  position  uncertainty 
on  timing  solutions  has  been  examined.  Results  based  on  the  assumption 
that  timing  errors  at  one  station  are  perfectly  known  indicate  the 
following: 

(1)  The  error  in  the  solution  of  time  bias  terms  is  correlated  witi 
radiation  pressure  modeling  errors  and  appears  to  be  approximately  a 
linear  function  of  its  magnitude. 

(2>  Multiple  satellite  processing  yields  a better  determination  of 
system  timing  parameters  and  satellite  orbital  parameters  than  single 
satellite  processing  in  the  presence  of  unmodelled  radiation  error.  In 
this  case  radiation  pressure  was  not  a parameter  in  the  solution. 

(3)  Adequate  timing  solutions  from  single  satellite  processing  are 
achievable  if  accurate  radiation  pressure  modeling  is  available. 

(4)  The  uncertainty  in  station  position  affects  the  standard  errors 
on  timing  solutions.  Therefore  station  position  should  be  known  to 
better  than  one  meter. 
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APPENDIX  A 
PASS  MATRIX  CONCEPT 
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PASS  MATRIX  CONCEPT 


Data  Aggr elation 

The  pass  matrix  concept  as  used  in  orbit  determination  involves 
aggregating  satellite  tracking  data  on  a pass  basis.  Normal  equations 
are  formed  from  the  tracking  data  of  each  separate  pass  over  a station. 
The  parameters  of  fit  are  the  satellite  orbital  constants  and  certain 
bias  parameters  characteristic  to  the  pass.  These  bias  parameters  re- 
present three  components  of  station  position,  refraction  bias,  range  or 
frequency  bias,  and  range  of  frequency  drift. 

Formation  of  Pass  Normal  Equations 

After  data  from  a pass  has  been  filtered  to  eliminate  "bad"  points, 
the  normal  equations  are  formed: 

Let  Dt.  represent  the  data  taken  at  time  t^  with  associated  standard 
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where  the  parameters  of  fit  Pj,...P^,  C^,  and  bias  terms  are  given  in 
Table  A-l.  The  weight  matrix  W for  the  data  of  a pass  is  given  by 


A-l 


1 


W3 


0 


(A- 2) 


With  the  A and  W matrices  defined,  the  least  squares  normal  equations 
for  a pass  are  given  by: 


BAP  = E 

(A- 3) 

B = ATWA 

(A-4) 

E = A^W6 

(A-5) 

where  the  n x 1 vector  6 contains  the  observational  residuals.  The  pair 
’Ej,  E " are  denoted  as  the  pass  matrices  for  the  j'th  pass. 
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PARAMETERS  OF  FIT 


NOTATION 


P4 


D 

Bias  Terms 
x,y,z 

CR 

RI!  °r  fB 
*B  °r 


PARAMETER 
x or  a 
y or  e sin  a 

z or  e cos  a. 

• 

x or  I 

y or  L+G 

» -\ 

z or 

Drag 

Station  Position 
Refraction  Bias 
Range  or  Frequency  Bias 
Range  or  Frequency  Drift 


TABLE  A-l 
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APPENDIX  B 
LEAST  SQUARES  METHOD 


B 


LEAST  SQUARES  METHOD  FORMULATION 


Range  Data 

The  range  data  class  is  given  by 

tut)  = | Rgat(t)  - T?g  (t)  | + Rg  + Rb  - t0>  + -Rd  «•  cR)  (B-n 
where  the  vectors  Rsat  and  Rs  are  the  radius  victors  of  the  satellite 
and  observing  station  respectively.  The  terms  R^  and  are  respectively 
the  range  bias  and  range  drift  characteristic  to  the  pass.  represent? 
a refraction  scaling  parameter  and  AR  the  Hopfield  tropospheric  refrac- 
tion correction.  The  time  tQ  is  the  epoch  of  the  fit. 

Least  Squares  Method 

Let  denote  the  pass  matrices  (Appendix  A)  formed  from  the 

filtered  observations  from  the  i*th  pass  of  satellite  j over  station  k. 
These  normal  equations 

= *1  (B'2^ 

are  formed  for  the  orbital  (7)  and  bias  (6)  parameters  associated  with 
the  satellite-station  pair  jk  (Table  A-l). 

Kith  each  timing  offset  modeled  as  a bias  and  linear  drift,  the 
range  and  range  drift  pass  parameters  for  the  i'th  pass  are  written, 
respectively,  as 


Rb  = (tj  + tk)  c 


n»-3> 


and 


R^  s (tj  + tfc)  C <R-4> 

where  tj  and  t^  are  the  time  bias  for  satellite  j and  station  k and  tj 
and  t^  are  the  time  drifts.  These  time  bias  parameters  are  not  however 


B-l 


pass  dependent. 

With  the  partial  derivatives  of  the  range  data  D(t)  with  respect  to 
time  bias  and  time  drift: 


ol)(t)  oD( t)  e c 


ot, 


(B-5) 


Mti  = c(t.to> 

Otk 


(B-6) 


**j 


the  pass  matrix  pair  be  expanded  with  the  four  time  bias 


parameters  replacing 


Bl’  Eljk“ > BhEi'jk. 

(13  parameters)  (15  parameters) 


(B-7) 


Kor  each  pass  i the  elements  of  the  pass  matrices  are  re- 

1 1 jk 

arranged  so  that  pass  dependent  parameters  (station  bias  and  refraction 
correction)  may  be  formally  eliminated: 
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(B-8) 


jk 


where  the  subscript  "b"  refers  to  pass  dependent  parameters  and  the 
subscript,  "o"  refers  to  all  other  parameters  (orbit,  drag  and  timing). 
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The  equations  for  the  elimination  of  pass  i dependent 

parameters  are 

Boo  <AVjk  + Bib  <AV  jk  = I’ 

(B-9a) 

Bbo  <“0>Jk  +Bib  ^“b)jk  = ^- 

(B-9b) 

Solving  equation  (B-9b)  for  (4Pb)jk  and  substituting  equation  (B-9a)  we 

obtain 

^Boo  “ BobB6b®bo)  ^po>jk  = Bo  " BobBbbBb 

(B- 10) 

or 

Bi  (i"Vjk  = Ei 

(B— 11) 

For  each  pair  jk  we  now  have  the  eliminated  pass  matrices  "bJ,  e”* 

for  the  parameter  set  CP  > fi.e  union  of  all  parameter  sets,  (??  ) 

jk  * ° jK 

consists  of  6j  orbital  parameters,  j drags,  and  2(j  + k)  timing 

parameters.  By  keeping  track  of  those  parameters  in  each  matrix  pair 

V ^i*jk’  the  elcments  fr,,m  al1  eliminated  pass  matrices  may  now  be 
combined  to  form  the  normal  equations  for  all  orbital  and  timing 
parameters. 
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